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functions.  It  is  assumed  that  the  generalized  response  is  a  nonlinear  function  of 
tbe  present  input  and  a  material  property-history  kernel  determined  by  a  Steklov- 
Lebesgue  norm.  Examples  applied  to  solid  propellant  are  given.  , 
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SUMMARY 

v 

A  broad  class  of  materials  possessing  both  instantaneous  nonlinear 
elasticity  and  dissipation  in  addition  to  fading  memory  with  aging  effectt. 
is  described.  The  measure  of  the  generalized  input  function.  A,  which 
is  a  multl-plet  in  F,  the  deformation  gradient;  0,  the  temperature; 
g  =  grade,  as  well  as  various  chemical  affinities,  A^;  is  given  by  a 
semi-norm  over  a  Banach  Space.  With  the  definition  of  the  history 
At  *  At(s)  s  (t-s);  se[0,®)  and  with  the  restriction  of  At  to  past 
history  given  by  A«*  =  Aj*(s}  b  Ajj(t-s);  se(0,«)  the  sani-norm  is: 
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(0)^(0) 


with 


1+k  measff 
id) 


■,|>p(s)dsj^p<  «  and  hp(s)  >  0  . 


and  AV1;(0)  the  present  1ch  time  rate  of  change. 


(1) 


The  second  summand  is  a  modification  of  the  Steklov  Average 
[Kantorovich,  1964]  to  P  -  integrable  Lebesgue  functions. 


It  is  assumed  that  the  generalized  response  ft(t)  is  a  nonlinear 
functional )  of  the  present  input  A(t)  and  a  material  property- hi  story 
kernel  determined  by  the  Steklov-Lebeegue  norm  UAHs)]  UN  such  that 

P 


(2) 


Equation  (1)  shows  that  as  time  Increases  mesD  increases  and 
that  the  Influence  of  the  past  history  on  present  response  decreases. 

For  a  given  finite  duration  inputs  ite  influence  deoreaeeo  both  the 
longer  in  the  past  it  occurred  and  the  older  the  material  is.  This 
latter  effect  is  termed  Steklcv-aging.  As  the  age  of  the  material 
becomes  very  large  the  past  history  effects  are  obliterated  and  from 
Eq.  (1) 

11m|  |A*(s)|  |SN  =  3>Sa^(o)  jh^io)  •  (3) 

£  -+■  oo  "P 

mesO  ■+■ «  1Bo 
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Thus  a.  Steklov-aging  material  under  lemg  term  aging  approaches 
the  behavior  of  a  nonlinear  viscoelastic  or  Markoffian  type  material. 
For  very  small  Inputs,  A(o)  and  small  time  rates  of  change  of  Inputs, 
AU/to),  the  material,  after  long  term  aging,  becomes  linearly 
viscoelastic. 

Examples  of  the  theory  as  applied  to  solid  propellants  and  a 
sand-asphalt  concrete  are  given. 


1-0  INTRODUCTION 


The  application  of  the  linear  and  rvonlinssr  theories  of  viscoelas¬ 
ticity  has  had  wide  application  in  the  characterization  and  stress 
analysis  of  solid  propellants,  highly  solids  loaded  polymers,  and 
asphalt  pavement  during  the  past  two  decades. 

So  long  as  the  characterization  and  verification  of  the  selected 
constitutive  equation  was  confined  to  the  usual  loading  or  straining 
conditions  such  as  a  constant  strain  rats,  ramp  relaxation,  or  constant 
stress  creep  the  comparison  of  theory  with  experiment  was  generally 
satisfactory.  This  result  was  to  be  expected  since  the  theory  and 
experiment  are  often  a  curve  fitting  exercise  over  a  narrow  range  of 
load  types. 

When  viscoelastic  theory  is  applied  to  the  prediction  of  repeated 
loads  such  as  a  sawtooth  input  or  multiple  ramp  input,  the  predictive 
results  of  the  theory  are  often  quite  unsatisfactory.  The  above  state¬ 
ment  applies  to  linear  viscoelasticity  as  well  as  to  the  several  forms 
of  nonlinear  viscoelasticity  utilizing  the  multiple  integral  approach. 

Details  of  the  previous  statements  are  presented  in  [Farris  and 
Fitzgerald,  1970],  [Fitzgerald  and  Farris,  1970],  ar,d  Chapter  XI  of 
[Fitzgerald  and  Hufferri,  1971].  A  most  excellent  comparative  review 
is  given  in  [Stafford,  1969]. 

Figure  1  demonstrates  the  above  problem  based  upon  some  of  Farris' 
experiments,  using  a  filled  polyurethane  propellant. 

The  use  of  a  repeated  sawtooth  -train  history  is  typified  by  the 
curves  of  Figure  2  for  a  highly  solids  loaded  polybutadiene  acrllo- 
nitrile  propellant  [Bennett,  1971]. 
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Previous  publications  by  the  author  and  his  co-workers  snployed 
the  use  of  homogeneous  functions  of  the  strain  history  made  specific 
through  Lebesgue  norms.  The  Lebesgue  norm,  jjejj 
weighted  P-summable  Integral  of  the  infinitesimal  strain  history 
wherein 


,  is  essentially  the 
LP 


/ 

o'rru 


|ejP'  hp(s)ds 
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meas  S 


The  term  hp(s)  Is  a  positive  decreasing  function  of  s  for  fading 
memory  theories  where  P  is  taken  as  unity.  Certain  boundedness  restric¬ 
tions  apply  to  hp(s). 

For  P  •*■ «,  the  Chebyshev  norm  results  wherein 


v 


|  |ej  ji  =  ess.  sup  |e] 

Lco 


(5) 


Using  a  rather  simple  expression  [Farris,  1970]  for  stress,  T, 
versus  strain  e,  wherein 


T11(t)  *  435 


11*1111. 12-25 
TT^TfiTj  . 


vn 


(t) 


(6) 


the  results  of  Figure  3  were  obtained  [Farris,  1970,]. 

It  should  be  noted  that  the  general  form  of  E a.  (6),  which  Is  the 
!.m  norm  divided  by  the  Lp  norm  to  then**1  power  (with  P  >  1),  when 
multiplied  by  e.jj  yields  the  following  results: 


from  which  P+1,  hence  P  and  n  may  be  determined. 

The  predictions  of  Figure  3  were  made  using  data  obtained  as  described. 
A  comparison  of  the  above  norm  expression  with  linear  viscoelasticity  for 
a  polyurethane  propellant  Is  given  In  Figure 4 [Farris,  1970]. 

A  more  general  expansion  of  Farris’  earlier  expressions  has  been 
given  by  [Vakily  and  Fitzgerald,  1972]  wherein  the  stress  function 
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involves  a  sum  of  Lebesgue  norm  ratios  and  the  naturally  occurring 
inverse,  power  law  kernel  in  the  viscoelastic  integral  as  follows: 


Tn(t)'=Z  L Afj ( TFlTp-)  10 

i=o  d«o  \  pi  /  0 


(t~T)"m3  e]](c)dr 


The  first  term  of  the  above  expansion,  i  s  j  =  0  with  q  =  m 


yields 


Tit \t)  =  A 


l«ii(t)IL> 


oo  e 


(t) 


which  is  the  previous  Farris  expression,  Eq.  (6), 

Including  the  current  value  of  the  strain,  |e-j-|j  and  simplifying 
a  three  term  expansion  of  Eq.  (10)  results  in 


Til (*)  =  A1 


i  I  inc 

hii 

w; 


£n(t)  +  ■  TFirrd 


(t-r)  e,i(r)dT  . 


where  we  have  defined  | |e| j0  a  je] . 

Various  other  specific  forms  of  the  above  Lebesgue  norms  may  be 


formulated. 


Applying  Eq.  (12)  to  a  series  of  compression  tests  on  a  sand- 
asphalt  mixture  produced  the  following  expression  [Vakily  and  Fitzgerald, 
1972]  where  the  constants  were  evaluated  from  constant  strain  rate 
"step"  relaxation  tests: 


The  moduli  for  a  strain  of  0.37%  and  0.5%  Is  shown  In  Figs.  5 
and  6. 

The  constant  strain  rate  tests  for  two  different  rates  are  given 
in  Fig.  /, 

The  ramp  relaxation  test  results  are  given  in  Fig.  8. 

Using  the  expression,  Eq.  (13),  derived  from  the  above  data, 
predictions  shown  as  open  circles  and  experiments,  shown  as  solid  lines, 
were  conducted  for 

(1)  an  interrupted  ramjj  strain  input  at  two  different 
strain  rates,  shown  in  Figs.  9  and  10.  A  comparison 
of  the  above  predictions  with  those  of  linear  visco¬ 
elasticity  is  given  in  Fig.  11 

(2)  a  single  constant  strain  rate  "sawtooth"  strain, 
shown  in  Fig.  12 

(3)  a  repeated  constant  strain  rate  test  cycled  between 
a  set  maximum  strain  and  zero  load,  Fig.  13. 

Again,  it  is  clear  that  the  sand-asphalt  material  exhibits  a  degree 
of  permanent  memory  similar  to  the  filled  polyurethane  propellant  and 
the  P8AN  propellant.  Figures  5  through  13  are  from  [Vakily  and  Fitz¬ 
gerald,  1972]. 

An  expression  for  the  stress  such  as 


T(t)  ■  f[U(t),  I|U<t>jjL  ] 

00 


(14) 
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produces  a  stress-strain  relation  with  no  relaxation  but  with  the  Mullins 
effect  predominant  [Mullins,  1947]. 

A  mathematical  justification  for  the  above  norms,  and  their  restric¬ 
tions,  is  given  In  [Hufferd  and  Fitzgerald,  1972]  including  thermodynamic 
implications. 

In  general,  one  may  express .the  Cauchy  stress  tensor,  T,  as  a 
function  of  the  norms,  for  example 

T-fOKt),  Hull,  ,-  Hull ,]  (15) 

—00  *-p 

where  U  is  the  positive  square  root  of  the  right  Cauchy-Green  (finite 
strain)  tensor.  The  resul  *?**-,•  polynomial  function  following  the  well 
known  Rivl In-Spencer  expv<.„>or  for  initially  isotropic  materials  will 
produce  an  expression  i.p  to  the  second  power  in  the  U  and  its  norms 
with  coefficients  that  are  polynomials,  or  preferably  here,  rational 
fractions  in  the  joint  invariants 

trU,  trU2,  -  -  ,  tr||U|j.  ,  etc. 

LP 

Again,. a  three-dimensional  expression  similar  to  Farris'  is  derived 
for  a  single  term,  homogeneous  form  for  T,  namely 

ftrl lull  f 

no  ” A  [trTnmTj  <ultl ' ”  (,6) 

Application  of  the  above  tensorial  equation  to  a  uniaxial  test 
produce  results  similar  to  those  previously  shown. 
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2.0  AGING  EFFECTS  AND  THE  STEKLOV  AVERAGE 


The  lebesgue  norms  previously  described  are  such  that  they  produce 
the  integral  of  the  strain  history  on  a  P  -  sumnable  basis.  The  use.  of  a 
non-unit  weighting  function,  hp(s),  will  provide  for  fading  .memory  effects. 
The  Lebesgue  norms  may  have  a  basis  in  microscopic  theory  as  presented  by 
Farris  in  his  doctoral  thesis. 

Nevertheless,  the  norms  as  used  herein  are  simply  continuum  postulates 
whose  use  is  justified  by  the  results  shown  herein  and  in  the  various  quoted 
references. 

With  the  integrals  as  used,  however,  ho  aging  effects  are  included. 

Consider  now  a  class  of  materials  whose  response  is  governed  by  certain 
weighted  averages  of  the  past  strain  history  as  well  as  by  the  present  value 
of  the  strain  and  its  several  time  derivatives. 

A  norm  on  such  a  space  may  be  constructed  as  follows  for  ?  generalized 
•input  A 


The  above  is  actually  a  semi -norm  unless  measD  is  finite, 

Considering  especially  the  second  summand  in  Eq.  (17),  we  shall  call 
it  a  steklov  Average  since  it  is  a  generalization  of  the  Steklov-Lebesgue 
average  gi ven i n  [Kantorovi ch  and  Akilov,  1964], 


A  slight  variation  in  the  above,  which  will  be  termed  the  modified  S;eklov 
Average  and  which  satisfies  all  the  requirements  of  a  norm  (or  semi -norm)  is 


with  k>0  . 


Dropping  the  fading  memory  factor  hp(s),  for  simplicity,  results  then  in 


where  we  have  taken  measD  =  t  for  the  usual  rather  smooth  physical  inputs,  A  . 

Looking  at  the  Lebesgue  ?}orm  ratios  of  Eq.  (10)  or  Eq.  (11)  for  example 
and  the  definitions  Eq.(4)  and  Eq,  (19)  yields  the  following  relation  between 
the  Lp  norms  and  the  $p  norms 


and  for  Q  -*•  »  , 
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It  is  readily  shown  that  the  ratio 

IMS 

WT?  -'  VQ>P'  P-’ 


since  the  dominance  of  Lq  over  Lp  is  well  known  when  Q>P  [Kantorovich  and 
Akilov,  1964]  and  fur  k>0,  the  multiplier  is  also  >1  . 

Defining  T1^  as  the  stress  obtained  from  Lebesguc  norms  as  In  the  previous 
section  1.0,  the  stress  obtained  by  substituting  modified  Steklov  norms  is 


T  where 


TS  =  TL  (1  +  kt)^~  n 


or  with  Q  -*■ 00 


TS  ■  TL(1  +  kt) 


Consider  for  example  a  material  governed  by  a  simple  permanent  memory 


norm  relation  such  as  Eq.  (6) 


I  I  nO 


L  .  ,p  1*11  II  1 
1  lji«iiiii. 


en  (0 


For  a  step  relaxation  test  with  strain  magnitude  z  (25)  yields 


t11L-  Aeot"n/P 


Applying  Instead  the  ratio  of  the  modified  Steklov  norms  Sp  produces 
Tii S  3  A  en  t"n/P  (1  +  kt*)n/P 

»  I  0 

,11. 


3 


(27) 


ic 

For  a  material  which  is  strained  immediately  upon  being  formed,  t  =  t  . 
Otherwise  however,  t  is  the  time  relative  to  the  beginning  of  the  step  strain 
whereas  t  is  the  actual  time  from  the  Initial  creation  of  the  material  to 
the  present.  The  difference  in  behavior  of  the  Lp  and  $p  norms  is  shown 
in  fig.  14.  A  step  load-unload  input  is  used  with  the  L-j,  S,  and  lro  =  = 

e  max.  curves  plotted  for  an  example.  After  the  load  occurs,  L-j  is  constant 
as  is  Lw.  However,  the  norm  reduces  with  time  since  it  is  averaging  over 
the  life  of  the  specimen. 

The  relation  between  the  Lp  and  Sp  norms  is,  from  (4)  and  (19) 

1  |A|  I-  s  1 1 A|  |  (1  +  kt)-,/P  (28) 

5P  LP  • 

It  will  also  be  noticed  from  Eq.  (24)-  that  the  exponent  of  the  "aging" 

if 

term  1+kt  is  equal  to  the  positive  numerical  slope  of  the  inverse  power 
law  relaxation  modulus,  n/P,. 

Consider  a  typical  solid  propellant  with  n/P  3  0.2  ,  then 

Tn  3  A  V”0,2  +  kt*)0,z  *  (29) 

If  the  relaxation  modulus  Increases  by  20%  over  a  one-year  ambient  aging 
then  k  =  3  x  10"8,  sec"^  . 

A  typical  CTPB  propellant  increases  its  modulus  by  20%  in  100  days. 

With  a  slope  n/P  3  0.2  ,  hence,  the  value  of  k  is  15  x  10  sec  . 

if 

It  is  also  to  be  noted  that  for  large  values  of  kt  ,  the  aging  is 
described  by  a  straight  line  on  a  log-log  scale.  Further  for  n/P  3  0,  no 

relaxation  and  no  strain  rate  effects,  the  aging  is  nonexistent. 
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Because  of  the  srnal  <  values  of  k,  the  use  of  the  Lp  rather  than  the  Sp 

norms  for  short  time  loads  cn  newly  formed  materials  Is  fully  equivalent. 

*  . 

Further,  for  short  loading  times  relative  to  the 'Tife  of  the  material, 
(l+kt*)n^P  is  essentially  a  constant  and  Steklov- aging  reduces  to  the  so- 
called  nomothetic  aging  process  [DeArriaga,  19691.  J 

For  finite  strain,  the  essential  results  also.  hold.  Consider  an  incom¬ 
pressible  material  subject  to  a  simple  step  elongation  .with  a  stretch  ratio 
of  X  where  X  is  the  stretched  length  divided  by  the, initial  length.  Using 

t  * 

a  simple  constitutive  equation  of  the  form  of  (16)  but  with  modified  Steklov 
norms  instead  of  Lebesgue  norms  as  in  (21)  results  in 

T(t)  =  t”n</P  (1  +  kt*)n/P  (X  -  1)  ‘  (30) 

For  a  material  stretched  when  newly  formed  and  then  held  at  the  stretch 
* 

ratio,  X,  with  t  =  t 


For  long  times-  then  as  t  +  «  , 

T(t)  «  k°/P  (x  -  1)  (32) 

Lx+2x~*J 

One  could  obviously  extend  the  expansion  to  higher  order  terms  In  U, 
buc  the  essential  point  to  be  made  is  that  the  material  does  not  fully  relax 
as  t  +  «  .  For  the  typical  material  previously  mentioned,  if  the  relaxation 
modulus  at  one  second  were  1000  psi,  the  final  modulus  would  reduce  to  60  psi . 

In  addition,  there  can  be  considered  an  elastic  component  with  no  loss  of 
generality.  ,13. 


3.0  CONCLUSIONS 


It  has  been  proposed  that  the  characterization  of  relaxing,  rate  sensitive 
materials  be  based  upon  certain  weighted  averages  of  the  past  history  which 
are  semi-norms  called  modified  Steklov  averages. 

The  rationale  is  based  upon  previously  mentioned  guidance  from  molecular 
theory  but  is  primarily  based  upon  the  postulate  that  the  present  response  of 
a  material  is  governed  by  Its  present  deformation  gradient  and  a  selected 
weighted  P-suiranable  Steklov  average  (17). 

Since  in  engineering  practice,  one  seldom  knows  the  detailed  past  of  a 
structures  thermal  and  deformation  history  but  usually  knows  the  average  and 
the  maxi mums,  no  serious  loss  of  applicability  should  result.  . 

The  justification  for  the  use  of  the  proposed  norms  is  based  upon  the 
several  examples  given  herein  where  the  predictions  using  norms  is  shown  to 
be  superior  in  accuracy  to  the  predictions  based  upon  linear  viscoelasticity. 
The  use  of  the  ^  or  maximum  value  norms  implies  that  the  material  also  is 
sensitive  to  the  maximum  strain  value  it  has  ever  been  subjected  to.  The 
possible  extension  of  the  concept  to  viscoplastic  materials  is  under  study. 

Farris  has  shown  that  for  a  selected  class  of  solid  propellants  the 
use  of  the  conventional  time- temperature  superposition  integral  is  valid. 

The  present  author  suggests  that  the  time  temperature  shift  integral 


for  redded  time  5  with  the  temperature  shift  factor  a T  be  used  in  the  modified 
Steklov  norm  for  aging  effects.  It  is,  of  course,  equally  possible  to  use  an 
absolute  reaction  rate  correction  for  the  kt  term  in  (19). 
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The  use  of  a  rather  simple  homogeneous  fraction  involving  only  the  max. 
norm  and  the  P-norm  of  the  modified  Steklov  average  (25)  produces  three 
specific  results  for  infinitesmal  deformations 

.the  relaxation  modulus  Is  described  by  an  Inverse  power  law,  t“n/,p 
.the  constant  strain  rate  curves  are  described  by  a  power  law  whose 
exponent  is  unity  plus  the  relaxation  modulus  exponent,  t1_n^P 
.the  relaxation  and  secant  moduli  are  subject  to  an  age  hardening 
process  n/P  described  by  the  factor  (1+kt) 

.for  very  long  term  stretching  and  aging,  the  material  stress  response 
is  nonvanishing  as  shown  In  (32). 

It  has  been  also  shown  that  for  short  aging  times  the  form  of  the 
equations  reduces  to  ratios  of  Lebesgue  norms  only.  For  short  term  loading 
relative  to  the  age  of  the  material,  a  homothetic  aging  process  is  produced. 

There  does  not  exist  a  unique  Inverse  for  strain  as  a  function  of  stress, 
i.e. ,  a  general  creep  Inverse.  Indeed,  an  Inverse  does  not  generally  exist 
unless  a  nonconstant  weighting  functiont  hp(s)  is  used.  Even  then,  there 
results  a  nonlinear  Integral  equation  of  the  type 

Xsf  -fmi  hp(s)ds  ;  *(5)  ■  AP(s)  (34) 

with  the  exponent  m  a  rational  fraction  whose  value  will  generally  be  near  unity. 

Since  the  present  value  of  stress  Is  governed  by  certain  averages  of 
the  past  strain  history,  uniqueness  in  the  Inverse  is  not  to  be  expected. 

It  Is  therefore  suggested  that  a  creep  law  be  formulated  In  the  same  fashion 
as  (25),  for  example 

e(t)  -  f[||T||j  .  I|T||S  .  T(t)]  (55) 

p  00 


“15” 


made  specific  as- 


e(t)  =  B 


|  IT  j  j  V 


rrm 


m 


T(t) 


which  for  a  constant  stress,  TQ,  results  in 


e(t)  =  B  T0  tm/P  (1  +  kt)"m/P 


(36) 


(37) 


The  calculation  time  necessary  to  use  the  various  norm  equations  given 
herein  is  much  shorter  for  general  inputs  than  the  time  needed  for  visco¬ 
elasticity.  This  computer  time  saving  results  from  the  fact  that  the  norms 
are  only  a  number  at  the  present  time  whose  value  changes  by  the  modified 
average  at  each  step.  If,  however,  nonconstant  weighting  functions,  hp(s), 
are  used,,  no  computational  advantage  results. 
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Figure  3  -  COMPARISON  OF  CALCULATED  AND  OBSERVED-  STRESS-STRAIN  OUTPUT  FOR  AN  INTERRUPTED  RAMP  STRAIN  INPUT  . 


Figure  4  -  STRESS  OUTPUT  FOR  INTERRUPTED  CONSTANT  STRAIN  RATE  TE 


Figure  5  -  RELAXATION  MODULUS  OF  SAND-ASPHALT  FOR  FIRST  TEST 


Figure  6  -  RELAXATION  MODULUS  OF  SAND-ASPHALT  FOR  SECOND  TEST 


STRESS  OUTPUT  FOR  TWO  DIFFERENT  CONSTANT  STRAIN  RATE  TESTS 


Figure  8  -  STRESS  OUTPUT  FOR  TWO  DIFFERENT  RAMP  RELAXATION  TESTS 
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Fioure  9  -  STRESS  OUTPUT  FOR  INTERRUPTED  RAMP  STRAIN  INPUT 


Calculated  Points  for  Non-linear  Viscoelasticity 


Figure  11  -  COf-PARISGN  OF  CALCULATED  AND  OBSERVED  STRESS-STRAIN 
OUTPUT  FOR  AN  INTERRUPTED  RAMP  STRAIN  INPUT 
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Figure  12  -  COMPARISON  OF  CALCULATED  AND  OBSERVED  STRESS-STRAIN 
BEHAVIOR  OF  SAND-ASPHALT 
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